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Abstract. We consider finite systems of interacting Brownian particles including 
active friction in the framework of nonlinear dynamics and statistical/stochastic 
theory. First we study the statistical properties for 1 — d systems of masses connected 
by Toda springs which are imbedded into a heat bath. Including negative friction 
we find A*' + 1 attractors of motion including an attractor describing dissipative 
solitons. Noise leads to transition between the deterministic attractors. In the case 
of two-dynamical motion of interacting particles angular momenta are generated 
and left/right rotations of pairs and swarms are found. 



1 Introduction 

The main purpose of this work is the study of the dynamics of active Brow- 
nian particles inchiding interactions. The friction is modelled by a velocity- 
dependent function derived from a model of energy supply The 
interaction between the particles is modelled by Toda potentials (in 1 — d) or 
Morse potentials (in 2 — d). 

Since the classical work of Toda the development of the nonlinear dynam- 
ics and statistical thermodynamics of of Toda systems has remained a central 
topic of research. Toda was able to find exact solutions for a special 1 — d sys- 
tem with a exponential interaction. In particular Toda proved the existence of 
soliton solutions and calculated the exact partition function. Since that soli- 
tons as excitations of nonlinear chains of masses found a remarkable interest. 
Several interesting results were obtained for the statistical thermodynamics 
of Toda systems We consider here systems with active friction and 

noise by means of analytical tools and simulations. It was shown in previous 
work, that there exists a special temperature regime around a transition tem- 
perature Ttr between the phonon and the soliton regime, in another context 
called the localization temperature Tioc, where the interaction of solitons and 
phonons is strong and has a remarkable influence on several physical phenom- 



ena, a special one being energy localization at special sites |13 and another 
one the excitation of a broadband coloured noise spectrum with an 1 // region 
at low frequencies iQ. Here our main interest is devoted to the influence of 
negative friction on the properties of Brownian motion. Driving the system 
by negative friction to far from equilibrium states we find + 1 attractors 
of deterministic motion including an attractor describing dissipative solitons. 
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Noise leads to transition between the deterministic attractors. In the case 
of two-dynamical motion of interacting particles positive or negative angular 
momenta are generated with equal probability. This leads to left/right rota- 
tions of pairs, clusters and swarms. We will show that the collective motion 
of large clusters of driven Brownian particles reminds very much the typical 
modes of parallel motions in swarms of living entities. 

2 Equations of motion, friction and forces 

Let us consider a systems of point masses m with the numbers 1,2, ...N. 
We assume that the mass m is connected to the next neighbours at both sides 
by Toda forces, The distance between the mass i and the mass i + 1 is de- 
noted by the equilibrium distance is assumed to be ct, therefore the spring 
elongation reads = Ri — a. In the following take a as the length unit. 

The dynamics of the system is given by the following equation of motion 
for the elongations 

^r, = [V'in+i) - 2V'in) + y'(r,_i)] - jiv,)v, + T,{t) (1) 

where Ti{t) is a stochastic force with strength D and a (5-correlated time 
dependence: 

{T^t)) =0; {T,{t)T,{t')) = 2DS{t - t')8,^ (2) 

In the case of thermal equilibrium systems we have 7(w) = 7(0) = const.. 
In the general case where the friction is velocity dependent we will assume 
that the friction is monotonically increasing with the velocity and that the 
limit for large velocities is a well defined constant: 

7(w) — > 7(0) = const. (3) 

Our basic assumption is, that in this limit the loss of energy resulting from 
friction, and the gain of energy resulting from the stochastic force, are com- 
pensated in the average. From this postulat follows the non-equilibrium fluct- 
uation-dissipation theorem : 

D = kBTjn/m (4) 

Here T is the temperature of the heat bath, ks is the Boltzmann con- 
stant, and D is the strength of the stochastic force. For the case of a passive 
thermal heat bath 7 = 70 eq.(^) agrees with the conventional Einstein rela- 
tion. The validity of a fluctuation-dissipation relation between noise strength 
and damping strength is assumed in order to guarantee the excistence of a 
stationary or thermal equilibrium, independent on the limit of the friction pa- 
rameter 7o > 0. We note that 70 = 0, T = corresponds to the conservative 
case. 
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In the following wc will study first the simplest case of passive friction, i.e. 
we will assume that the friction function is constant "f(v) = 70 = const. Then 
more complicated friction functions will be studied. Historically velocity- 
dependent friction forces were first studied by Rayleigh and Helmholtz. Ex- 
tensions of these models were investigated in many papers on driven Brow- 
nian dynamics [|l5j. A characteristic property of these friction functions is 
the existence of a zero of friction for a finite velocity vq which defines a kind 
of attractor in the velocity-space. The model we use here for active friction 
with an attractor for the velocities vq was introduced in 0] and j|] . Detailed 
studies of this so-called energy depot model model may be found in |^,||. 
The case of Rayleigh- friction was analyzed in detail in [Q . The depot-model 
of the friction function is based on a concrete model of Brownian motion with 
energy supply, storage in a tank and conversion into motion 0] . We assumed 
that the Brownian particle itself is capable of taking up external energy stor- 
ing some of this additional energy into an internal energy depot, e{t). Within 
the energy depot model the depot may be may be altered by three different 
processes: 

1. take-up of energy from the environment; where q is the pump rate of 
energy 

2. internal dissipation, which is assumed to be proportional to the internal 
energy. The rate of energy loss is denoted by c. 

3. conversion of internal energy into motion, where d is the rate of conversion 
of internal to kinetic degrees of freedom. The depot energy is used to 
accelerate the motion. 

Our model of energy supply is motivated by investigations of active bi- 
ological motion, which relies on the supply of energy. The supplied energy 
is in part dissipated by metabolic processes, but can be also converted into 
kinetic energy. The energy depot model leads in an adiabatic approximation 
to the following friction function 

^(^) - ^0 - YTd^ 

For the energy depot model of active friction exists an attracting velocity. 
This is defined by the zero of the friction function 7(^0) — which has the 
value 

yl^^;a^^^-l (6) 

Here a plays the role of the bifurcation parameter. For a < the friction is 
purely passive, i.e. in average no energy is supplied. For a > the friction 
function has a negative part near to u = and a zero point which acts as an 
attracting set in the velocity space. 

The interaction between the masses i and i -|- 1 is in the Id-case modelled 
by the Toda force 

/, = -V^'(rO = y(exp[-H-l) (7) 
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Here b is the stiffness of the springs and u is the hnear oscihation frequency 
around the equiUbrium position. 

In the harmonic hmit the force is given by fi — ujri . The spring energy is 
described by the Toda potential 

9 

TflUJ 

V{n-uo,h) = ^{ex:p[-hri\^l + hn) (8) 
We will use mcj^cr^ a,s the energy unit. 

In general the calculation of the distribution functions of interacting parti- 
cles is not a trivial task since it is connected with the solution of a multi- 
dimensional Fokker-Planck equation. A full solution is available only for the 
passive case 7 = 70 = const. Then the only attractor of the dynamics is 
the rest state of the particles and the statistical properties in equilibrium are 
described by a canonical Toda ensemble under pressure P and temperature 



T |10| , |liyi2| . The stationary solution of the Fokker-Planck equation which 
)ur Langevin equation reac 

PaiPn, rn) = exp ( -- 



corresponds to our Langevin equation reads 

p2 /2m + Veff 



2—1 ^ (9) 

^ V27rmfci3r exp {X) r{X + Y) ^ ' 

with the effective potential 

Veff{r,,)^V{rn)+Prn (10) 

Here fcs denotes the Boltzmann constant, and V{r) is the Toda potential 
defined above. Further we used the abbreviations 

The elementary excitations in passive Toda rings including the noise spec- 
trum and the structure factor were investigated in detail in the work . 
The investigation of the dynamics of Toda rings with energy supply by ac- 
tive friction is more difhcult and only partial solutions are available |,|. We 
discuss here the special case when the supply of energy is given by the de- 
pot model in the approximation of a velocity-dependent friction [0J^,^ . Due 
to the driving slow particles are accelerated and fast particles are damped. 
At definite conditions our active friction function has a zero corresponding 
to stationary velocities wq, where the friction disappears. The deterministic 
trajectory of our system is then attracted by a cylinder in the 4d-space given 

by 

vl + vl - vl (12) 



where Vq is the value of the stationary velocity which is for the depot-model 

a 
d ■ 



given by Vg 
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In Fig. [I] we see the result of a histogram calculated by J. Dunkel for N = 
8 particles with Toda interactions Here the result of 1000 runs with 
stochastic initial conditions is represented as a function of the average velocity 
of the particles. The shape of the histogram demonstrates the existence of 
A'' + 1 attractors of motion. We underline however that the histogram shown 
in Fig. 1 is the result of finite time runs, so it represents the shape of a finite- 
time distribution function. In the limit of long runs the distribution may 
change; we expect e.g. that at long times the statistical weights are shifted 
from the center of the distribution to the wings. 

3 Two-dimensional dynamics 

The effects demonstrated in the previous sections are not restricted to 1 — d 
Toda lattices but persist at least qualitatively in more realistic 2 — d and 3 — c? 
models of dense fluids consisting of solvent and solute molecules with Morse- 
or Lennard- Jones interactions. Superposition of solitons corresponds to mul- 
tiple collisions in these systems. In higher dimensions a weak localisation of 
potential energy was observed also at the bindings of the bath molecules 
and was connected to a transition between different lattice configurations 

We introduce interactions described by the potential U{ri, ...,rAr); then the 
dynamics of Brownian particles is determined by the Langevin equation: 

vecr^=Vt; mVi ^ -j{vi)vi - WUiri, ...,rN) + ^{t) (13) 

where J^(U) is a stochastic force with strength D and a (5-correlated time 
dependence as defined above. 

We will discuss now the motion of active particles in a two-dimensional space, 
r = {xi,X2}. The case of constant external forces was alredy treated by 
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Fig. 1. Probability distribution for 8 active Brownian particles 
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Schienbein and Grulcr pc| ] . Symmetric parabolic external forces were studied 
in 1^,^,^ and the non-symmetric case is beeing investigated in 1 2l| . 
Here we will study 2d-systems of > 2 particles (see Fig. |2 ) 
We will begin with the study of a 2 particle problem . Let us imagine 
two Brownian particles which are pairwise bound by a pair potential U{ri — 
r2)- The two molecules will form dumb-bell-like configurations. Then the 
motion consists of two independent parts: The free motion of the center of 
mass having the coordinates Xi = 0.5(a;ii -|- X21) and X2 = 0.5(xi2 -|- X22)- 
The relative motion under the influence of the potential is described by the 
coordinates xi = (xn — X12) and X2 = {xi2 — X2i). The motion of the center 
of mass M is described by the equations 



Xi = Vi MVi = -7 (^1,^2) Vi + Ti{t) 
±2 = V2 MV2 = -7 {Vi,V2) V2 + T2{t) 



(14) 



The stationary solutions of the corresponding Fokker-Planck equation 
reads [B| 



P^{V)^C{l + dV^)^'''^''^ exp 



2D 



(15) 



This corresponds to the driven motion of a free particle located in the 
center of mass. 

The relative motion is described by the equations 



ii — vi fj,vi = —7 (wi, U2) vi — diU + Ti{t) 

X2^V2 fJ,V2 = ~^ ivi,V2) V2 - d2U + J^2it) 



(16) 




Fig. 2. Snapshot of 1000 Brownian particles rotating in a parabolic self-consistent 
field 
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with /i = m/2. For simplification we specify now tlic potential U{r) as a 
symmetric parabolic potential 

U{xi,X2)^^a{xl + xl) (17) 

First, we restrict the discussion to a deterministic relative motion, which 
is described by four coupled first-order differential equations. The relative 
motion of 2 particles corresponds to the absolute motion of 1 particle in an 
external field. Therefore we will omit now for simplicity the tilde denoting the 
relative character of the motion and denoting the mass again by m instead 
of ^ : 



Xl — Vi mil — —7 (^^1, "^2) "^1 — axi 

±2=^2 mv2 = {VI,V2)V2 ~ ax2 (18) 

For the one-dimensional case it is well known that this system pocesses 



a limit cycle corresponding to sustained oscillations 1 15 . For the 2 — d case 
we have shown in Q that a limit cycle in the 4 — d space is developed, which 
corresponds to left/right rotations with the frequency uq. The projection of 
this periodic motion to the {vi, V2} plane is the circle 

vj +vl ^ vl ^ const. (19) 
The projection to the {xi,X2} plane also corresponds to a circle 

xl+xl=rl^^ (20) 

The energy for motions on the limit cycle is 

Eo^j{vl+vl) + ^{xl+xl) (21) 

- Y^o + (22) 

We have shown in Q , that any initial value of the energy converges (at 
least in the limit of strong driving) to 

H — > Eq = mvl (23) 

This corresponds to an equal distribution between kinetic and potential 
energy. As for the harmonic oscillator in 1 — d, both parts contribute the 
same amount to the full energy. This result was obtained in based on 
the assumption that the energy is a slow (adiabatic) variable which allows a 
phase average with respect to the phases of the rotation. In explicite form 
we may represent the motion on the limit cycle in the A ~ d space by the 4 
equations 
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xi = ro sm{LUot) ; vi = — tq ujq cos{ujQt) (24) 

2^2 = fQ cos(woi) ; ''^2 = — foujQ sm{ujt) (25) 

The frequency is given by the time the particle need for one period moving 
on the circle with radius tq with constant speed vq. This leads to the relation 

1/2 



(26) 



ujo 

This means, the particle oscillates (at least in our approximation) with 
the frequency given by the linear oscillator frequency uj. 

The trajectory on the limit cycle defined by the above 4 equations is like 
a hula hoop in the 4 — d space. The projections to the xi — X2-space as well 
as the projections to the vi — i;2-space are circles. The projections to the 
subspaces X1-V2 and X2-V1 are like a rod. In the 4 — d space the attractor 
has therefore the form of a hula hoop. A second limit cycle is obtained by 
reversal of the velocity. This second limit cycle forms also a hula hoop which is 
different from the first one, however both I.e. have the same projections to the 
{xi,X2} and to the {vi,V2} plane. The motion in the {xi,X2} plane has the 
opposite sense of rotation in comparision with the first limit cycle. Therefore 
both limit cycles correspond to opposite angular momenta. L3 = +/LtroWo and 
L3 = —firaVQ. Applying similar arguments to the stochastic problem we find 
that the two hoop-rings are converted into a distribution looking like two 
embracing hoops with finite size, which for strong noise converts into two 
embracing tyres in the A — d space. In order to get the explicite form of the 
distribution we may introduce the amplitude-phase representation |^ 

xi = psin{ujQt + 4>)\ vi ~ ~ pujQCOs(ujQt + (p) (27) 

X2 — pcos{ujot + 0) ; V2 ~ — pujQsin{ujQt + 0) (28) 

where the radius p is now a slow and the phase </> is a fast stochastic 
variable. By using the standard procedure of averaging with respect to the 
fast phases we get the distribution of the radii . 



Po(p) - C (1 + dpu,^)'-'"''^ exp [-^pc^o^J (29) 

The probability crater is located above the two deterministic limit cycles 
on the sphere Tq = Vq/ujq. Strictly speaking not the whole sherical set is filled 
with probability but only two circle-shaped subsets on it, which correspond 
to a narrow region around the limit sets. The full stationary probability has 
the form of two hula hoop distributions in the 4d space. This was confirmed 
by simulations 

The projections of the distribution to the {xi,X2} plane and to the 
{vi, V2} plane are smoothed 2(i-rings. The distributions intersect perpendicu- 
lar the {xi,V2} plane and the {x2,vi\ plane. Due to the noise the Brownian 



70 2 



Active Brownian Particles 9 



particles may switch between the two hmit cycles, this means inversion of the 
angular momentum (direction of rotation) [p|j21 



Summarizing our findings for a 2-particle system forming a dumb-bell 
system we may state: The center of mass of the dumb-bell will make a driven 
Brownian motion corresponding to a free motion of the center of mass. In 
addition the dumb-bell is driven to rotate around the center of mass. What 
we observe then is a system of rotating Brownian molecules. The internal 
degrees of freedom are excited and we observe driven rotations. 

In this way we have shown that the mechanisms described here may be 
used also to excite the internal degrees of freedom of Brownian molecules. 

An extension of the theory of pairs leads to a theory of the motion of 
clusters of active molecules 1 22 2^ , ^ , p4[ . In Fig. ^ we see a snapshot of sim- 
ulations of the stochastic dynamics of a cluster of 1000 active interacting 
Brownian particles. 

In order to simplify the simulations we assumed that the interaction of 
the molecules in the cluster is given by a van der Waals - like interaction with 
a relatively long range range. For example we may use the interaction model 
proposed by Morse 



(r) = A [exp(-ar) - if - A 



(30) 



Due to the attracting tail the molecules form clusters. The individual molecules 
move then in the coUectice (self-consistent) field of the other molecules which 
might be represented by a mean field approximation 



V{r) = / dr'(t){f - r')p{r') 



(31) 



where r = (ii,i2) is the radius vector counted from the center of mass and 
p(r') is a mean density in the cluster. Approximating V by quadratic terms 




Fig. 3. Snapshot of the same swarm after one half period of rotation 
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only we get 

V{xuX2) = "^0 + ^ {aiij + 02^1) + ... (32) 

In this way we arrive again at the harmonic problem we have studied 
above. We may conclued that the individual molecules in the cluster move 
at least in certain approximation in a parabolic potential. From this follows 
that they will make rotations in the field. We have performed simulations 
with 1000 particles moving in a self-consistent potential of parabolic shape. 
The driving function has a zero at Wq = 1. In Figs. ^ and |^ one sees two 
snapshots of the moving cluster formed by the molecules. 

Since the indidual particles move in an effective parabolic potential an- 
gular momenta are generated, the swarm starts to rotate. Similar as in the 
case of the dumb-bells the clusters are driven to make spontaneous rotations. 
Finally a stationary state will be reached which corresponds to a rotating 
cluster with nearly constant angular momentum (see Fig. U). Under the in- 
fluence of noise the cluster may switch to the opposite angular momentum, 
i.e. to the opposite sense of rotation, the system occurs to be bistable. 

In Fig. 4 is shown that in noisy systems of active Brownian particles 
the two values L3 = —mrovo, +mroVo have the maximal probability. Strong 
coupling of the particles leads to synchronization of the angular momenta, 
for weak coupling the cluster may be decomposed into groups with different 
angular momentum |2l||2^ . The rotating swarms simulated in our numerical 
experiments remind very much the dynamics of swarms studied in papers of 



Viscek and collaborators |22,E3| and in other recent work fel 24 



4 Conclusions 

We studied here the active Brownian dynamics of a finite number of particles 
including interactions modelled by Toda or Morse potentials and velocity- 
dependent friction. Beside analytical investigations we have made a numer- 
ical study of relatively small systems. For this we have been investigat- 
ing Id-model systems with nonlinear Toda interactions and 2(i-models with 
parabolic confinement. At low enough temperatures we observe only har- 
monic excitations and the models reduce to systems of coupled harmonic os- 
cillators. However at higher temperatures the dynamics is completely changed 
and nonlinear excitations of soliton-type come into play. At high tempera- 
tures the behaviour of Id interacting systems is dominated by soliton-type 
excitations. Including active friction we observe in the Id-case -t- 1 attrac- 
tors of nonlinear excitations. For interacting Brownian particles in 2d-case 
angular momenta are generated (see Fig. 4) and we find rotations of the par- 
ticles around the center of attraction and collective rotational excitations of 
swarms. 
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Fig. 4. Probability distribution of the angular momentum for active Brownian 

tides with Vq = Vq = 1 



